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A new method of treating electronic energy level transitions as well as linking ionization 
to electronic energy levels is proposed following the particle-based chemistry model of Bird. 
Although the use of electronic energy levels and ionization reactions in DSMC are not new 
ideas, the current method of selecting what level to transition to, how to reproduce transition 
rates, and the linking of the electronic energy levels to ionization are, to the author’s 
knowledge, novel concepts. The resulting equilibrium temperatures are shown to remain 
constant, and the electronic energy level distributions are shown to reproduce the Boltzmann 
distribution. The electronic energy level transition rates and ionization rates due to electron 
impacts are shown to reproduce theoretical and measured rates. The rates due to heavy 
particle impacts, while not as favorable as the electron impact rates, compare favorably to 
values from the literature. Thus, these new extensions to the particle-based chemistry model 
of Bird provide an accurate method for predicting electronic energy level transition and 
ionization rates in gases. 
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Nomenclature 

= diameter of particle (m) 

= number of transitions in sample 
= sample time (s) 

= distribution function 
= energy level degeneracy 
= energy levels 

= Boltzmann’s Constant (k= 1.3806503e-23 m 2 kg s -2 K" 1 ) 
= forward reaction rate (m 3 /molecule/s) 

= orbital quantum number 

= total orbital angular momentum quantum number 
= mass (kg) 

= number of molecules 

= number density (1/m 3 ) and principal quantum number 
= probability 

= total spin quantum number 
= temperature (K) 

= relative velocity in a collision (m/s) 

= number of active electrons for an atomic level 
= Dirac delta 
= energy (J) 

= translational energy ( Str = 0.5 p v r 2 ) 

= gamma function 

= reduced mass in a collision ( p = mimflmHmi ) ) 

= exponent in variable hard sphere model 
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0 

= characteristic temperature (K) 


Subscripts 

0,1 

= ground and first electronic energy level 

coll 

= collision value 

diss 

= dissociation value 

el 

= electronic energy component 

heavy 

= heavy particle (atom/molecule) 

i,j,J 

= energy levels 

ion 

= ionization value 

max 

= maximum value 

pail- 

= value obtained from both particles in collision 

particle 

= value obtained from single particle in collision 

ref 

= reference value 

tr 

= translational energy component 

vib 

= vibrational energy component 


I. Introduction 

As we consider missions, both manned and unmanned, where the (re)entry velocities result in shock-layer 
temperatures on the order of tens-of-thousands of degrees Kelvin, the importance of electronic energy levels and 
ionization reactions become more pronounced. Although the treatment of electronic energy levels and ionization 
reactions using the direct simulation Monte Carlo (DSMC) method is not a new development 1-3 , most of these 
methods are based on measured, equilibrium rates, which are always questionable at high temperature, especially 
when applied to non-equilibrium problems. The primary reason for the inadequate state of chemical reaction 
modeling is the difficulty in accurately measuring the internal energy state specific reaction rates to validate 
theoretical models in the temperature range of interest. Therefore, most reaction rates are based on low-temperature 
equilibrium measurements and are fit to the most reliable measured data sets. In some cases, the uncertainty 
associated with the spread of these measurements exceeds one order of magnitude 4 . 

Recently, an approach for determining chemical reaction rates from microscopic molecular data, referred to as 
the particle-based chemistry model (PBM), has been developed 5, 6 that does not use any macroscopic rate 
information. It is the purpose of this paper to introduce a new method of treating electronic energy level transitions 
following the PBM methodologies. One distinct feature of the PBM is the direct linking of the vibrational energy 
level of a molecule to the dissociation of that molecule. This idea is extended in this paper to the direct linking of 
the electronic energy level of an atom/molecule to ionization. 


II. Detailed Electronic Energy Level Model for DSMC 

In order to implement a model for the distribution of electronic energy in the DSMC technique, there are three 
procedures that must be defined. First, when a particle is first introduced into a simulation, it is necessary to obtain 
a new electronic energy from a given distribution through statistical sampling. This is also required when a particle 
changes its electronic energy after a collision with a surface. Under equilibrium conditions, the distribution has the 
well known Boltzmann form. Second, it is necessary to statistically sample a new electronic energy following a 
collision which involves electronic energy transfer. Finally, a method to reproduce the electronic energy transition 
rate must be defined. In the following subsections, each of these tasks are considered. 

A. Equilibrium Sampling 

Each electronic energy level j has a distinct energy, Sj, and degeneracy, gj. The Boltzmann distribution for the 
electronic energy levels at a given temperature T gives the following result for the fraction of particles in level j: 


( 1 ) 


N, 8)exp(-e,-/g) 
N 

£g,exp(-e,./^r) 

1=0 


This distribution is used when creating a new particle in the DSMC simulation at a boundary specified at the 
temperature T. However, it is not possible to sample an electronic energy level j directly from the distribution. 
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Therefore, an acceptance-rejection procedure is used. This is performed by selecting values for j from the following 
distribution: 


( 2 ) 


f , ( g,eM~ e J lkT ) 

gj exp(-£ y / kT ) 


where J is the value of j for which Eq. ( 1 ) is a maximum. Unfortunately, since the degeneracy of each level is a 
variable specific to each species, there is no way to know a priori what level is going to give the maximum, so the 
maximum level must be searched for at each implementation or saved for a constant boundary temperature. The 
sampling of a new electronic energy level then proceeds as follows: 

1. select at random an electronic energy level evenly distributed between 0 and Jmax, where Jmax is the 
maximum possible energy level. A random number, RANDi, evenly distributed between 0 and 1 is needed such 
that: 

(3) t = int[(l + J mjt ) * RAND, ] 

2. determine the value J for which Eq. (1) is a maximum; 

3 . accept the value of j if f e i ’(j) > RAND 2 ; 

4. if the value of j is not accepted, then return to step 1. 


B. Post-collision sampling 

A phenomenological approach is usually adopted in the DSMC method when a collision occurs that involves 
energy transfer. The Borgnakke-Larsen method 7 samples a post-collision state from a combined distribution of the 
translational and electronic collision energies of the colliding particles. Based on the approach of Bergemann and 
Boyd 8 and Boyd 9 , the Dirac delta is used to write the distribution of energies in Eq. (1) in the following continuous 
form: 


(4) 


fe,{ £ elj) 


gj™ P(- £ j /kT ) 

l ma\ 

X^ ex P(~ e / 1 kT ) 



1=0 
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Consideration must now be given to the distribution of translational energy of the colliding particles. This 
distribution is naturally affected by the intermolecular model used since this determines the collision probability. 
For the present study, the variable hard sphere (VHS) collision model of Bird 10 is used. Elowever, it is a simple 
matter to develop the formulation for an alternative collision model such as the variable soft sphere of Koura 11 . For 
the VHS model, the distribution of the translational collision energies is: 


(5) 



1 

r(y 2 -co) 



exp 



Using Eqs. (4) and (5), the combined distribution for sampling a post-collision electronic energy level j ' from the 
total collision energy e co ii = £tr + s e i = s V + s ’ e i is: 


(6) 




r ( 5 A- (0 ) g i ex p(~ e / / kT ) 
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In applying the general Borgnakke-Larsen scheme, it is assumed that local thermodynamic equilibrium prevails. 
Therefore, the temperature T in Eq. (6) is constant. Also, the total collision energy is constant, so it is only 
necessary to perform sampling of the post-collision state from the following distribution form: 


(7) 


8(f\Ecou) x 8y(l 


'coll 



Again, an acceptance-rejection procedure is used. The normalized distribution that is required is obtained by finding 
the value of j’ for which Eq. (7) is maximum, J\ which is different for each value of e C oii. Therefore, the following 
distribution is obtained: 


( 8 ) 
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In the DSMC code, these procedures are implemented in the following way: 

1. Given a pair of particles with total collision energy e co ii which undergoes electronic energy exchange, 
determine J’; 

2. determine the maximum allowable electronic energy level obtainable from e C oiu J”', 

3. take J * to be the smaller of J’ and J”\ 

4. as described in the previous section on equilibrium sampling, perform an acceptance-rejection procedure to 
sample j ’ from: 


(9) 
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C. Electronic Energy Level Transitions 

A transition from electronic energy level i to level j for atomic nitrogen, N, can be written in the form of a 
chemical reaction as: 

(10) N 1 + M — > N j + M 

Following the PBM, the simplest model is to assume that the transition occurs if the electronic energy level of N 
after a trial Borgnakke-Larsen redistribution of the collision energy is j* =j. Since we are treating the transition as a 
chemical reaction, there are Jmax transitions to consider where Jmax is the maximum energy level energetically 
possible. There is a possibility that there will be multiple transitions where j * = /, each of which, according to the 
PBM, are equally probable. In the DSMC code, the decision to transition from level i to j is implemented in the 
following way: 

1. check each transition i to j from j = 0 to Jmax, performing a trial Borgnakke-Larsen redistribution for each, 
keeping a list of the possible reactions (j * =j). 

2. if there are any possible transitions from step (1), choose a random transition from the list and perform the 
transition. 


III. Linking Ionization to Electronic Energy Levels in DSMC 

Bird 4 5 proposed a set of molecular-level chemistry models based solely on fundamental properties of the two 
colliding molecules: their total collision energy, quantized vibrational energy levels, and molecular dissociation 

energies. These models link chemical-reaction cross sections to the energy exchange process and the probability of 
transition between vibrational energy states. The Borgnakke-Larsen procedures and the principle of microscopic 
reversibility are then used to derive simple models for recombination reactions and for the reverse (exothermic) 
reactions. These models do not require any macroscopic data, and they function by seeking to balance the fluxes 
into and out of each state, thus satisfying microscopic reversibility. 
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Table 1. Gas properties. 


Test Gas 

cl ref (m) 

m (kg) 

T ref (K) 

ra 

©vib (K) 

Zref 

0diss (K) 

£ion (J) 

N2 

3.580E-10 

4.650E-26 

1000 

0.68 

3371.0 

52500 

113500 

2.496 192E- 18 

N 

3.110E-10 

2.325E-26 

1000 

0.65 

- 

- 

- 

2.327964E-18 

e 

5.640E-15 

9.109E-31 

1000 

0.60 

- 

- 

- 

- 


In this section, these ideas are extended to include the quantized electronic energy levels and molecular/atomic 
ionization energies. For dissociation to occur, if the collision energy: 


( 11 ) 


’'coll 


F + F 

c 'tr,pair ' c 


vib , particle 


exceeds the dissociation energy, the molecule in question dissociates. If this is extended to electronic energy levels 
and ionization, if the collision energy now defined as: 


( 12 ) 


"'coll 


— F + F 

c 'tr,pair ' ° 


el , particle 


exceeds the ionization energy, the particle in question loses an electron, becoming an ion. 

IV. Results 

In this section, the new electronic energy exchange, electronic energy level transition, and ionization reaction 
models are applied in a two-dimensional DSMC code where the results are accumulated over all cells to 
approximate a zero-dimensional simulation. The test gas is comprised of molecular and atomic nitrogen, Nj and N, 
and electrons, e, where the relevant simulation parameters are listed in Table 1. A list of electronic energy levels for 
N 2 and N are presented in Tables A. 1 and A.2 in the Appendix, respectively. The simulations were run at varying 
temperatures between 10,000 K and 60,000 K at a number density of le23/m 3 in an adiabatic box 0.002 m on a side 
with 300,000 molecules. During the simulations, relaxation of the particles was allowed to proceed as usual, but 
when a reaction was determined to take place, the number of reactions was advanced by one but the simulators were 
left unchanged, so no energy was added to or taken away from the flow. In the figures included in this section, the 
solid lines are values from quoted rates from the literature, and the symbols are the values that were sampled from 
DSMC. 


A. Electronic Energy Level Transitions and Comparison to Equilibrium 

In order for these procedures to be acceptable, they need to reproduce equilibrium conditions after the simulation 
has been allowed to run for a sufficient amount of time. This set of 
simulations were allowed to proceed for approximately 24 hours, 
resulting in at least 300 millions collisions. The first requirement is 
for the procedures to remain at the input equilibrium temperature. 

The input equilibrium temperature is compared to the sampled 
translational temperature in Fig. 1 and can be seen to remain 
constant. The resultant population distribution over the first twenty 
energy levels for a range of temperatures is compared to the 
Bolzmann distribution in Fig. 2. Again, the simulation represents 
the equilibrium state very well. 

Sampled transition rates will now be compared to those found in 
literature. We must first define how to measure the sampled 
transition rates in DSMC. The transition rate for the reaction shown 
in Eq. (10) is calculated as: 


dn . 


(13) 


k f 


dt 

n N > n N 



Sampled T (K) 

Figure 1. Sampled translational 
temperature over a range of input 
equilibrium temperatures. 
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Figure 2. Sampled electronic energy level distributions for the first 20 levels of atomic nitrogen over a range 
of equilibrium temperatures. 


where the numerator is the change in the number density of level i due to the reaction in the sampled time. 

1. Electron impact rates 

We will begin by examining electronic energy level transitions caused by electron impacts. Sampled transition 
rates are compared to those by Chemyi 12 in Fig. 3 for molecular nitrogen. The sampled rates very nearly match the 
values quoted in the literature. Sampled electron impact electronic energy level transition rates for atomic nitrogen 
are compared in Fig. 4 with values from Frost 13 with similar results. 



- N° -> N 3 

• 

N° -> N 3 (DSMC) 

- W°->W J 

■ 

N° -> N 4 (DSMC) 

- N 1 -> N 3 

♦ 

N 1 -> N 3 (DSMC) 

- A/’-s/V 4 

▲ 

N 1 -> N 4 (DSMC) 

- N 2 -> N 3 

T 

N 2 -> N 3 (DSMC) 
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Figure 3. Electron impact electronic energy 
transition rates for N 2 . 
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Figure 4. Electron impact electronic energy 
transition rates for N. 


6 

American Institute of Aeronautics and Astronautics 






, 





0 10000 20000 30000 40000 50000 6000C 

Temperature (K) 

Figure 5. Electron impact ionization rates for N 2 . 
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Figure 6. Electron impact ionization rates for N. 


Next, electron impact ionization rates will be compared for both atomic and molecular nitrogen. Figure 5 
presents a comparison between the sampled rates and those from Losev 14 for molecular nitrogen while Figure 6 
compares the rates between DSMC and a variety of sources 4, 14-19 for electron impact ionization rates for atomic 
nitrogen. As with the electronic energy level transition rates for electron impact, the ionization rates compare well 
with literature. 

2. Heavy particle impact rates 

Using the current model for electronic energy level transitions and ionization, the comparison between rates 
measured using DSMC and those from the literature for heavy particle impacts are not as similar as those for 
electron impacts. Sampled heavy particle impact endothermic electronic energy level transition rates for molecular 
nitrogen are compared to those of Park 20 in Fig. 7. The measured trends match those from the literature for 
individual level transitions, but the magnitudes are different, especially for the level 0 to 1 transition for a Ni + N 2 
collision. The other transition rates are approximately within an order of magnitude of the quoted rates, but the 
comparison is not as good as the electron impact transition rates. Sampled heavy particle impact exothermic 
electronic energy level transition rates for N 2 are compared to those of Guerra 21 in Fig. 8. The quenching transitions 
are much closer to the measured rates as compared to the excitation transitions. This is most likely due to the fact 
that it is the quenching rates that are measured in the laboratory. The excitation rates are then calculated assuming 
equilibrium. It is possible that the transitions in the laboratory are not at equilibrium, therefore some of the rates are 
farther off than others. 

When examining atomic nitrogen heavy particle impact electronic energy level transition rates, only quenching 
transition rates were found for comparison to the sampled rates. The sampled rates are compared to values quoted 
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Figure 7. Heavy particle impact electronic energy 
level endothermic transition rates for N 2 . 
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Figure 8. Heavy particle impact electronic energy 
level exothermic transition rates for N 2 . 


1 

American Institute of Aeronautics and Astronautics 


by Levin 22 in Fig. 9. Once again the quenching rates for 
electronic energy transitions compare reasonably well. Heavy 
particle impact ionization rates are not prevalent in the 
literature, so are not included in this analysis. 


V. Conclusion 

Kinetic theory based chemical reaction models recently 
proposed for the DSMC method have been extended to include 
electronic energy level transitions as well as ionization. This 
new set of models do not use measured macroscopic reaction/ 
transition rates to calibrate adjustable parameters. Instead, they 
make use of the principles of microscopic reversibility and 
molecular-level energy exchange to predict the probability that 
a chemical reaction or energy level transition occurs during a 
collision between two particles. 

Procedures have been defined that are required in order to 
implement the current model for the distribution of electronic 
energy in the DSMC technique. These include defining 
methods to sample from an equilibrium (Boltzmann) 
distribution, to sample post-collision states from a combined 
distribution of the translational and electronic collision 
energies of the colliding particles, and to determine when a particle should transition electronic energy levels. The 
first two methods are straight forward and can be derived following methods for other energy modes. The last 
method is a novel implementation of the new chemical reaction models. These three methods have been shown to 
maintain equilibrium temperatures and reproduce the Boltzmann distribution of electronic energy levels after 
sufficient simulation time. Sampled electronic energy transition rates have been shown to compare well to values 
from literature for electron impact transitions and reasonably well for heavy particle impact transitions. 

The new chemical reaction models have also been extended to include ionization reactions. Following the new 
dissociation model, where a molecule dissociates if the available collision energy exceeds the dissociation energy, a 
particle will loose an electron if the available collision energy exceeds the ionization energy. The sampled 
ionization rates have been shown to agree well with values from literature for electron impact reactions. However, 
heavy particle ionization data are scarce in literature and have been excluded from the current study. 

Appendix 

Table A.l. Listing of electronic energy level information of atomic nitrogen. 


Level j 

r'w,/ (J) 

Si 

0 

0.00000 

1 

1 

9.9726950E-19 

3 

2 

1 . 1 843030E-1 8 

6 

3 

1 . 1 880540E- 1 8 

6 

4 

1.3164670E-18 

3 

5 

1 .3538 160E-1 8 

1 

6 

1.3762700E-18 

2 

7 

1.4482750E-18 

2 

8 

1.5414820E-18 

5 

9 

1.6924570E-18 

1 

10 

1.7242390E-18 

6 

11 

1.7706570E-18 

6 

12 

1.8473930E-18 

10 

13 

1.9388470E-18 

6 

14 

2.0778240E-18 

6 


8 



Temperature (K) 

Figure 9. Heavy particle impact 
electronic energy level exothermic 
transition rates for N. 
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Table A.2. Listing of electronic energy level information of atomic nitrogen. 


Level j 

Belj (J) 

gj 

»j 

k 

Xj 

Si 

Li 

Config. 

Term 

0 

0.00000 

4 

2 

1 

3 

3/2 

0 

2s 2 2p 3 

4 S° 

1 

3.8195294E-19 

10 

2 

1 

3 

1/2 

2 

2s 2 2p 3 

2 D° 

2 

5.7287486E-19 

6 

2 

1 

3 

1/2 

1 

2s 2 2p 3 

2po 

3 

1.65541 72E- 1 8 

12 

3 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)3s 

4 P 

4 

1.7121 737E-1 8 

6 

3 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)3s 

2 P 

5 

1.7507040E-18 

12 

2 

1 

4 

3/2 

1 

2s2p 4 

4 P 

6 

1.8589476E-18 

2 

3 

1 

1 

1/2 

0 

2s 2 2p 2 ( 3 P)3p 

2 S° 

7 

1 .883901 9E-1 8 

20 

3 

1 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)3p 

4 D° 

8 

1.8972523E-18 

12 

3 

1 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)3p 

4po 

9 

1.9219038E-18 

4 

3 

1 

1 

3/2 

0 

2s 2 2p 2 ( 3 P)3p 

4 S° 

10 

1.9235458E-18 

10 

3 

1 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)3p 

2 D° 

11 

1.9426247E-18 

6 

3 

1 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)3p 

2po 

12 

1.9797645E-18 

10 

3 

0 

1 

1/2 

2 

2s 2 2p 2 (*D)3s 

2 D 

13 

2.0598225E-18 

12 

4 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)4s 

4 P 

14 

2.0697973E-18 

6 

4 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)4s 

2 P 

15 

2.0783847E-18 

6 

3 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)3d 

2 P 

16 

2.0801974E-18 

28 

3 

2 

1 

3/2 

3 

2s 2 2p 2 ( 3 P)3d 

4 F 

17 

2.0827249E-18 

12 

3 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)3d 

4 P 

18 

2.0828065E-18 

14 

3 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)3d 

2 F 

19 

2.0859128E-18 

20 

3 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)3d 

4 D 

20 

2.0884332E-18 

10 

3 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)3d 

2 D 

21 

2.115 1235E-1 8 

2 

4 

1 

1 

1/2 

0 

2s 2 2p 2 ( 3 P)4p 

2 S° 

22 

2. 1219883E-1 8 

20 

4 

1 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)4p 

4 D° 

23 

2.1257740E-18 

12 

4 

1 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)4p 

4po 

24 

2.1299658E-18 

10 

4 

1 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)4p 

2 D° 

25 

2.1343488E-18 

4 

4 

1 

1 

3/2 

0 

2s 2 2p 2 ( 3 P)4p 

4 S° 

26 

2. 1377136E-1 8 

6 

4 

1 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)4p 

2po 

27 

2.1827900E-18 

12 

5 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)5s 

4 P 

28 

2. 1 86701 7E-1 8 

6 

5 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)5s 

2 P 

29 

2.1896284E-18 

6 

4 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)4d 

2 P 

30 

2. 1 90 1 0 14E-1 8 

28 

4 

2 

1 

3/2 

3 

2s 2 2p 2 ( 3 P)4d 

4 F 

31 

2.191 64 1 9E- 1 8 

12 

4 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)4d 

4 P 

32 

2.191 64 1 9E- 1 8 

14 

4 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)4d 

2 F 

33 

2.1920256E-18 

14 

4 

4 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)4f D 

2[ 3 ]° 

34 

2.1927445E-18 

14 

4 

4 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)4f G 

2 [3]° 

35 

2.1930233E-18 

20 

4 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)4d 

4 D 

36 

2.1930782E-18 

18 

4 

4 

1 

1/2 

4 

2s 2 2p 2 ( 3 P)4f G 

2[4]° 

37 

2. 193 1259E-1 8 

10 

4 

4 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)4f D 

2 [2 ]° 

38 

2.1942236E-18 

6 

4 

4 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)4f D 

2 [1]° 

39 

2.1942468E-18 

10 

4 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)4d 

2 D 

40 

2.1944891E-18 

22 

4 

4 

1 

1/2 

5 

2s 2 2p 2 ( 3 P)4f G 

2 [5]° 

41 

2.1947443E-18 

10 

4 

4 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)4f F 

2[ 2 ]° 

42 

2.1949910E-18 

14 

4 

4 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)4f F 

2 [3]° 

43 

2.1950570E-18 

18 

4 

4 

1 

1/2 

4 

2s 2 2p 2 ( 3 P)4f F 

2 [4]° 

44 

2.1957241E-18 

10 

3 

1 

1 

1/2 

2 

2s 2 2p 2 (*D)3p 

2 D° 

45 

2.1992586E-18 

14 

3 

1 

1 

1/2 

3 

2s 2 2p 2 ( 1 D)3p 

2po 

46 

2.2061644E-18 

2 

5 

1 

1 

1/2 

0 

2s 2 2p 2 ( 3 P)5p 

2 S° 

47 

2.2090956E-18 

6 

5 

1 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)5p 

2po 

48 

2.2091778E-18 

20 

5 

1 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)5p 

4 D° 

49 

2.2092125E-18 

12 

5 

1 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)5p 

4po 

50 

2.2149140E-18 

4 

5 

1 

1 

3/2 

0 

2s 2 2p 2 ( 3 P)5p 

4 s° 

51 

2.2225265E-18 

10 

5 

1 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)5p 

2 D° 

52 

2.2310009E-18 

6 

3 

1 

1 

1/2 

1 

2s 2 2p 2 ( 1 D)3p 

2po 

53 

2.2374954E-18 

12 

6 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)6s 

4 P 
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Level j 

'f'V/,/ (J) 

gj 

»j 

l, 

Xj 

Si 


Contlg. 

Term 

54 

2.2391604E-18 

6 

6 

0 

1 

1/2 

l 

2s 2 2p 2 ( 3 P)6s 

2 P 

55 

2.2407743E-18 

6 

5 

2 

1 

1/2 

l 

2s 2 2p 2 ( 3 P)5d 

2 P 

56 

2.2408504E-18 

28 

5 

2 

1 

3/2 

3 

2s 2 2p 2 ( 3 P)5d 

4 F 

57 

2.2412458E-18 

14 

5 

3 

1 

1/2 

3 

2s 2 2p 2 ( 3 P 0 )5f 

2 [3]° 

58 

2.2414379E-18 

12 

5 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)5d 

4 P 

59 

2.2418528E-18 

14 

5 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)5d 

2 F 

60 

2.2420758E-18 

14 

5 

3 

1 

1/2 

3 

2s 2 2p 2 ( 3 Pi)5f 

2[3]° 

61 

2.2422583E-18 

18 

5 

3 

1 

1/2 

4 

2s 2 2p 2 ( 3 Pi)5f 

2[4]° 

62 

2.2423070E-18 

10 

5 

3 

1 

1/2 

2 

2s 2 2p 2 ( 3 Pi)5f 

2 [2 ]° 

63 

2.2428934E-18 

20 

5 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)5d 

4 D 

64 

2.2434580E-18 

10 

5 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)5d 

2 D 

65 

2.2436230E-18 

6 

5 

3 

1 

1/2 

1 

2s 2 2p 2 ( 3 P 2 )5f 

2 [ 1 ]° 

66 

2.2437593E-18 

22 

5 

3 

1 

1/2 

4.5 

2s 2 2p 2 ( 3 P 2 )5f 

2 [9/2]° 

67 

2.2438446E-18 

10 

5 

3 

1 

1/2 

2 

2s 2 2p 2 ( 3 P 2 )5f 

2 [2 ]° 

68 

2.2439954E-18 

14 

5 

3 

1 

1/2 

3 

2s 2 2p 2 ( 3 P 2 )5f 

2 [3]° 

69 

2.2440307E-18 

18 

5 

3 

1 

1/2 

4 

2s 2 2p 2 ( 3 P 2 )5f 

2 [4]° 

70 

2.2515070E-18 

14 

6 

1 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)6p 

4 D° 

71 

2.2522583E-18 

10 

6 

1 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)6p 

4po 

72 

2.2665626E-18 

12 

7 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)7s 

4 p 

73 

2.2679643E-18 

14 

6 

3 

1 

1/2 

3 

2s 2 2p 2 ( 3 P 0 )6f 

2 [ 3 ]° 

74 

2.2680191E-18 

6 

6 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P|6d 

2 P 

75 

2.2680197E-18 

6 

7 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)7s 

2 P 

76 

2.2682030E-18 

28 

6 

2 

1 

3/2 

3 

2s 2 2p 2 ( 3 P|6d 

4 F 

77 

2.2686007E-18 

12 

6 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)6d 

4 P 

78 

2.2688525E-18 

14 

6 

3 

1 

1/2 

3 

2s 2 2p 2 ( 3 Pi)6f 

2 [3]° 

79 

2.2689594E-18 

18 

6 

3 

1 

1/2 

4 

2s 2 2p 2 ( 3 Pi)6f 

2 [4]° 

80 

2.2689951E-18 

10 

6 

3 

1 

1/2 

2 

2s 2 2p 2 ( 3 Pi)6f 

2 [2]° 

81 

2.2690781E-18 

14 

6 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P|6d 

2 F 

82 

2.2699818E-18 

20 

6 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)6d 

4 D 

83 

2.2702472E-18 

10 

6 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P|6d 

2 D 

84 

2.2704316E-18 

6 

6 

3 

1 

1/2 

1 

2s 2 2p 2 ( 3 P 2 )6f 

2 [ 1 ]° 

85 

2.2705093E-18 

22 

6 

3 

1 

1/2 

5 

2s 2 2p 2 ( 3 P 2 )6f 

2 [5]° 

86 

2.2705473E-18 

10 

6 

3 

1 

1/2 

2 

2s 2 2p 2 ( 3 P 2 )6f 

2 [2]° 

87 

2.2706401E-18 

14 

6 

3 

1 

1/2 

3 

2s 2 2p 2 ( 3 P 2 )6f 

2 [3]° 

88 

2.2706620E-18 

18 

6 

3 

1 

1/2 

4 

2s 2 2p 2 ( 3 P 2 )6f 

2 [4]° 

89 

2.2753854E-18 

8 

7 

1 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)7p 

4 D° 

90 

2.2838695E-18 

12 

8 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)8s 

4 P 

91 

2.2843827E-18 

6 

7 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)7d 

2 P 

92 

2.2846102E-18 

6 

8 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)8s 

2 P 

93 

2.2846711E-18 

28 

7 

2 

1 

3/2 

3 

2s 2 2p 2 ( 3 P)7d 

4 F 

94 

2.2846925E-18 

10 

7 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)7d 

4 P 

95 

2.2854244E-18 

14 

7 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)7d 

2 F 

96 

2.2862873E-18 

20 

7 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)7d 

2 D 

97 

2.2864728E-18 

10 

7 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)7d 

2 D 

98 

2.2939594E-18 

6 

9 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)9s 

2 P 

99 

2.2949085E-18 

12 

9 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)9s 

4 P 

100 

2.2950 12 IE-1 8 

6 

8 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)8d 

2 P 

101 

2.2955 104E-1 8 

24 

8 

2 

1 

3/2 

3 

2s 2 2p 2 ( 3 P)8d 

4 F 

102 

2.296041 IE-18 

14 

8 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)8d 

2 F 

103 

2.2966880E-18 

12 

8 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)8d 

4 P 

104 

2.2967069E-18 

10 

8 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)8d 

2 D 

105 

2.2968674E-18 

20 

8 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)8d 

4 D 

106 

2.301 1420E-18 

6 

10 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)10s 

2 P 

107 

2.3020713E-18 

6 

9 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)9d 

2 P 

108 

2.3020392E-18 

20 

9 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)9d 

4 D 

109 

2.3023276E-18 

14 

9 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)9d 

2 F 
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Level j 

(J) 

gj 

«/ 


Xj 

Si 

Li 

Contlg. 

Term 

110 

2.3024486E-18 

12 

10 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)10s 

4 P 

111 

2.3039586E-18 

10 

9 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)9d 

2 D 

112 

2.3040740E-18 

12 

9 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)9d 

4 P 

113 

2.3063971E-18 

6 

11 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)lls 

2 P 

114 

2.3067336E-18 

12 

11 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P) 1 1 s 

4 P 

115 

2.3073584E-18 

6 

10 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)10d 

2 P 

116 

2.3074386E-18 

14 

10 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)10d 

2 F 

117 

2.3075347E-18 

20 

10 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)10d 

4 D 

118 

2.3090407E-18 

10 

10 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)10d 

2 D 

119 

2.3094253E-18 

12 

10 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)10d 

4 P 

120 

2.3098102E-18 

2 

3 

0 

1 

1/2 

0 

2s 2 2p 2 ( 1 S)3s 

2 S 

121 

2.3103385E-18 

6 

12 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)12s 

2 P 

122 

2.310466 7E-18 

12 

12 

0 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)12s 

4 P 

123 

2.3 1 125 1 7E-1 8 

6 

11 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)lld 

2 P 

124 

2.31 14120E-18 

14 

11 

2 

1 

1/2 

3 

2s 2 2p 2 ( 3 P)lld 

2 F 

125 

2.3 1 1 5722E-1 8 

20 

11 

2 

1 

3/2 

2 

2s 2 2p 2 ( 3 P)lld 

4 D 

126 

2.3129340E-18 

10 

11 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)lld 

2 D 

127 

2.3130302E-18 

12 

11 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)lld 

4 P 

128 

2.3135589E-18 

6 

13 

0 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)13s 

2 P 

129 

2.3142478E-18 

6 

12 

2 

1 

1/2 

1 

2s 2 2p 2 ( 3 P)12d 

2 P 

130 

2.3 158179E-18 

12 

12 

2 

1 

3/2 

1 

2s 2 2p 2 ( 3 P)12d 

4 P 

131 

2.3166991E-18 

10 

12 

2 

1 

1/2 

2 

2s 2 2p 2 ( 3 P)12d 

2 D 
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